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ABSTRACT 

We use numerical simulations of large scale structure formation to explore the cos- 
mological properties of Gamma-Ray Burst (GRB) host galaxies. Among the different 
sub-populations found in the simulations, we identify the host galaxies as the most 
efficient star-forming objects, i.e. galaxies with high specific star formation rates. We 
find that the host candidates are low-mass, young galaxies with low to moderate star 
formation rate. These properties are consistent with those observed in GRB hosts, 
most of which are sub-luminous, blue galaxies. Assuming that host candidates are 
galaxies with high star formation rates would have given conclusions inconsistent with 
the observations. The specific star formation rate, given a galaxy mass, is shown to 
increase as the redshift increases. The low mass of the putative hosts makes them 
difficult to detect with present day telescopes and the probability density function of 
the specific star formation rate is predicted to change depending on whether or not 
these galaxies are observed. 

Key words: cosmology: large-scale structure of the Universe - galaxies: formation - 
galaxies: evolution ~ gamma rays: bursts 



1 INTRODUCTION 

For three decades Gamma-Ray Bursts (GRBs) remained an 
enigma. Th e detection of optical afte rglows associated with 
a GRB (see Ivan Paradiis et al.ll200(l for a review), revolu- 
tionized our understanding of the phenomena. A GRB is 
now thought to occur when the core of a massive star col- 
lapses to form a black hole llWooslevl[T9 9^JPa£ZYngki|[l998l: 



iMacFadven fc Wooslevlll999l) fsee also iMeszarosI l200l for 

a review). In most cases when an optical afterglow from a 
GRB has been detected, they have been shown to occur 
in galaxies at intermediate or high redshifts. To date, the 
redshift of over 40 GRB afterglows and their host galax- 
ies have been measured. Because of the short lifetime of 
a massive star, a GRB is generated essentially instanta- 
neously after the star's formation, as compared to the evo- 
lutionary timescale of a typical galaxy or the cosmological 
timescale. The bursts are therefore generally considered to 
be good tracers of the his tory of massive star formation (e.g. 
iBlain Nataraianllioool) . 

It was early realized that because of their extreme 
brightness, GRBs and their afterglows might be detectable 
to very high redshifts and therefore might be useful as 
cosmological probes, in particular if th e bursts possessed 
a 'st andard candle' like property (e.g. iLamb fc Reicharg 
|200Q'1 . The main interest in the cosmological application 



of GRBs has resulted from the increased understanding 
of t he relation between GRBs and powe rful supernovae 
(e.g. iHiorth et aT]l2003l : IStanek et a'i]|2003^ and from there 
the mapping of the cosm ic star formation history (e.g. 
Ir a mirez-Ruiz et aTl |2002| . and references therein). Other 
possible applications have also been considered, such as the 
use of bursts as probes of the interstellar medium in evolv- 
ing galaxies and as probes of primordial star formation and 
reionization, although these are observation ally very chal- 
lenging at present (e.g. lDiorgovski et al .''2003'). Constructing 
a reliable Hubble diagram with GRB aft erglows, however, i s 
likely to remain a very difficult subject llBloom et alJl2003l) . 

In this paper, we take a different approach altogether. 
Instead of focusing on the diverse properties of the bursts 
themselves we use the fact that in almost all cases where 
sufficiently deep observations have been carried out, a host 
galaxy has been discovered. Many of the host galaxies are 
very faint, with R band magnitudes down to about 30 
Jjaunsen et a U Effll . and might have gone undetected had 
not a GRB occurred in them. We will occasionally refer 
to the hosts as GRB selected galaxies. Unfortunately, the 
number of hosts is still too limited for statistical stud- 
ies, but this is likely to change after Swift is launched 
in 2004. A number of first results on the hosts does ex- 
ist however, and a picture of their properties is beginning 
to emerge. The host morphology shows most of them to 
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be dwarf galaxies, either compact, irregular or interacting, 
possi bly merging, and they also tend t o be blue in color 
fe^Rloom ettjUToQa iT.e Floc'h et allEonl iRerger et alJ 
2003). Many hosts ar e sub-luminous, with L/L* < 1, (e.j 



Le Floc'h etai] 12003 With GRB 990705 a notable excep- 
tion (|Le Floc'h et al.ll2002h . The star formation rate {SFR) 



in most hosts as inferred from optical observations show low 
to intermediate activity, with SFR = 1 — bQMc^/yx (e.j 



iRloom et al.ll. 199d:ICa.stro-Tira.do et aLlEooHlFvubo et alJ 
I2OOII: ICharv et alJl2002l: iPrice et alJl2002h . In a number of 
cases the sub-mm and radio observations indicate a SFR 
that i s up to a factor of 1 higher than the optically inf erred 
rates jBe rgcr ct al."2003'). Finallv. ISokolov et al.l J200J) and 
ICharv et al., (.2002) infer a host galaxy mass, M, typically in 
the range 2-10*— 4-10^°Mq. In addition, the variation in the 
specific star formation rate, SFR/M, fr om host to host is 
modest, normally within a factor o f a few l|Oiaxy_etLa lJ20o3; 
IChristensen et al]l2004^ . In fact, IChristeiMer^^jTl20M) 
has shown that all the hosts in their sample have spectral 
energy distributions similar to young starburst galaxies with 
small to moderate extinction. Comparing their sample with 
the Hubble Deep Field, they find the hosts to be similar 
to the bluest of the field galaxies which have the highest 
specific SFR. 

As the bursts occur in galaxies, they can be viewed not 
only as tracers of massive star formation, but also as tracers 
of galaxy formation and evolution. The study of GRBs and 
their host galaxies is therefore intimately linked to the study 
of galaxy evolution in general. Although the presently known 
host sample indicates that bursts can occur in any type of 
galaxy, the majority of them are faint and show a dwarf-like 
morphology. They may therefore belong to the faint end 
of the galaxy luminosity function. GRR selected galaxies 
may prove to be an important sub-population of the galaxy 
population and their study may thus provide information 
about the galaxy luminosity function not easily accessible 
otherwise. 

The aim of this paper is to investigate theoretically the 
properties of the host galaxies of GRRs as well as their cos- 
mological evolution rather than to focus on the GRR pro- 
genitor and its properties. To this end, we use numerical sim- 
ulations of large scale structure formation with the aim of 
identifying objects in the simulations with properties corre- 
sponding to those observed in the h osts. The simulations we 
use, were previously considered by lAlimi fc Courtvl (l2004l) 
to explore a number of galaxy properties such as the cosmic 
star formation rate density and the galaxy mass function. 
The simulations reveal a significant population of low-mass 
galaxies whose existence raises a number of questions on 
galaxy evolution. In this paper we concentrate on the in- 
stantaneous and the specific star formation rates of the sim- 
ulated galaxies. We identify a sub-population of the galaxies 
with high specific star formation rate as the likely candidates 
for GRR hosts. We find that the majority of this putative 
host population has low-mass and low to modest SFR, al- 
though a number of high-mass objects also have similar spe- 
cific rates. We then explore the properties and cosmic evo- 
lution of the simulated host population and compare to the 
galaxy population as a whole. 

In section |21 we describe the details of the numerical 
simulations relevant for this study. In section|21we introduce 
the instantaneous star formation rate and in particular the 



specific star formation rate as a measure of the efficiency of 
star formation. We also discuss the properties of the general 
galaxy population at low redshift. In section 01 a similar dis- 
cussion is presented for high redshift and the cosmological 
evolution of the galaxy properties is explored. In section |S| 
we propose interpretations of the observations and we argue 
that GRR host galaxies may be a powerful tracer of the faint 
galaxy population. Section |S| concludes the paper. 



2 NUMERICAL SIMULATIONS 

The simulations were performed with a 3D N-body/hydro- 
dynamical code, coupling a PM scheme for computing grav- 
itational forces with an Eulerian approach for solving the 
hydrodynamical equations. Sh ock heating is treated with 
the a rtificial viscosity method jVon Neumann fc Richtmverl 
I1950I) . The radiative cooling processes included are: col- 
lisional excitation, coUisional ionization, recombination, 
bremsstrahlung and Compton scattering. CoUisional ion- 
ization equilibrium is not assumed and the cooling rates 
are computed from the evolution of a primordial compo- 
sition hydrogen-helium plasma. The simulations also in- 
clude a model for galaxy formation, summarized below. 
The numerical procedu res ar e described in more de tail in 
ICourtv fc AUmil J2004) and lAlimi fc Courtvl J2004) . The 
galaxy properties discussed in the latter paper were found 
to be consistent with observational data. This suggests that 
the galaxy formation model used here to investigate the na- 
ture of host galaxies of GRRs, captures the essence of the 
galajcy formation process. 

Numerical simulations of large scale structure formation 
do not currently allow following the formation of objects of 
mass lower than the scale of a grid cell, a few times 10^ Mq . 
But computing the thermodynamic properties of the gas and 
its cosmological evolution allow us to consider the physical 
conditions needed to form a galaxy. The most important 
condition is that the gas cloud is collapsing, in the sense 
that the cooling time is less than the dynamical time or 
the free fall time dRees fc Ostrikedri977l) . We then identify, 
at each time step, gas regions satisfying this condition. A 
fraction of the gas mass is turned into an object, labeled a 
"stellar particle" , with a rate: 



dniE 



'TT' 



(1) 



with niBit) the available baryonic mass enclosed within the 
gas region at cosmic time t, and t, a characteristic timescale. 
The cosmological evolution of these particles is followed in 
a non-coUisional way. To express the condition t^ooi < ta, 
we define the cooling timescale tcooi, computed from the in- 
ternal energy variation of the gas, and the dynamical time 
or free fall time, tg = 37r/32Gp. To make sure that gas 
regions giving birth to galaxies are correctly identified we 
use additional criteria: the size of the gas cloud must be less 
than the Jean's length given by Aj = Cs^Tv/Gp)^''^; the gas 
must be in a converging flow: V -v < 0; the baryonic density 
contrast, Sb = (5p/p)B, must be higher than a threshold 
(1 + 5b)s, here taken to be the value of the baryonic den- 
sity contrast at the turnaround, i.e. 5.5 eis co mputed in the 
top-hat collapse spherical model ijPadmanab han 1993). The 
total matter density, including dark matter, baryonic mat- 
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ter and stellar particles, is used in the expressions for the 
dynamical time and the Jean's length. In each cell, checking 
that the four criteria described above are satisfied, a frac- 
tion At/t, of the gas is turned into a stellar particle, with 
At the timestep of the simulation. If a stellar particle starts 
forming at cosmic time to when the baryonic mass enclosed 
within the grid cell is ms(to), we assume that the particle 
is fully formed at time ti = to + At. Using eq. we find 
that the final mass of the particle is given by 

rui = niBito) - ruBiti) ~ niBito)^- (2) 

The characteristic time is taken to be = max(tfi, 10* yr). 
The stellar particles are involved in the computation of the 
gravitational potential and their evolution is treated in the 
same way as the coUisionless dark matter. 

Galaxy-like objects are then defined, at any redshift, 
by grouping the stellar particles with a friend-of-friend algo- 
rithm. This algorithm joins together all particles closer than 
a distance proportional to a link parameter tj, usually taken 
to be 0.2 (equivalent to selecting an overdensity of 187 in 
the spherical collapse model). The mass range of the stellar 
particles spans from 10^ to 10* but in what follows 

we only consider, at any redshift, galaxies with a total mass 
higher than 5 • 10* Mq, although a number of lower mass 
galaxies also form in the simulations. At z = 0, this lower 
galaxy mass limit corresponds to a minimum of 11 stellar 
particles per galaxy and at z = 3 it corresponds to a min- 
imum of 27 particles. The total number of galaxy objects 
in the catalogs at z — 0, 1 and 3 are 1650, 2164 and 1602, 
respectively. 

Through the formation of the stellar particles and their 
groupings we therefore have access to the formation history 
of each galaxy at any redshift. The mass of the galaxy at any 
time is the sum of the mass of its individual stellar particles 
at that time. This mass is likely to be an estimate of all the 
stellar populations formed over the lifetime of the galaxy, 
as the dark matter component is not included in the galaxy 
mass. We estimate the epoch of formation for a given galaxy 
as the mean epoch of formation of all stellar particles the 
galaxy is composed of, weighted by the mass of each particle, 
i.e. tmi = '^(rriiti) / '^rrii, where ti is the cosmic time cor- 
responding to the epoch of formation of stellar particle i. An 
alternative method would be to define the formation epoch 
as the epoch, tm2, at which the galaxy has acquired half of 
its final mass. Both methods give similar epoch estimates as 
shown in Figure As the former method is based on the 
formation history of the stellar populations that make up 
the galaxy, it resembles the way galaxy ages are estimated 
observationally, namely from the age of the stellar popula- 
tions in the galaxy. In this work, all formation epochs are 
therefore estimated with the first method. 

The results of this paper are given for a A— cold dark 
matter model. The parameters of the simulations are: Ho = 
70 km s-i Mpc-\ Q,k = 0, fi™ = 0.3, = 0.7, 
fib = 0.02/i~^ with h = Ho/IQQ. The initial density fiuc- 
tuation spectrum uses the transfer functions taken from 
[Bardccn e t gd. J iggg.) with a shape parameter given by 
ISugivarnaT liliosl) . The fluctuation sp ectrum is normalized 
on COBE data (|Bimn fc Whitelll9^7^ leading to a fihered 
dispersion at _R = 8 Mpc of crs = 0.91. The resolution 
of the simulations is as follows: the number of dark mat- 
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Figure 1. Comparing the two different ways of estimating the 
epoch of formation of galaxies, tmi and tm2, computed from the 
galaxy catalog at z = 0. See text for definitions. 



ter particles is Np = 256"^ and the number of grid cells is 
Ng = 256^; the comoving size of the computational volume 
is Lio:c = 32 h-^Mpc, giving a dark matter particle mass 
of 2.01 X 10* Mq, and a gas mass initially enclosed in each 
grid ceU of 3.09 x lO'' Mq. 



3 INSTANTANEOUS STAR FORMATION 
RATE AND EFFICIENCY OF STAR 
FORMATION 

For each object at any given redshift z, corresponding to 
cosmic time t^, we compute the instantaneous star forma- 
tion rate SFR*, where * indicates a value obtained from 
the simulations. The instantaneous rate is computed from 
the mass in stellar populations in the galaxy that have age 
less than r at time t^, SFR* — AM/t. This is equivalent 
to the baryonic mass AM, having been converted to stars 
during time r prior to the epoch under consideration. Note 
that the SFR* depends both on r and the cosmic time, 
SFR* {T,tz). It depends also on the amount of gas in the 
galaxy that is available for star formation at tz since the 
star formation rate in the model is proportional to the bary- 
onic density. We use r = 0.1 Gyrs as representative value 
at z = 0, time dilated at higher redshifts. Figure |5| shows 
a comparison between SFR* computed with r — 0.05 and 
r = 0.1 Gyrs for galaxies in the catalog at 2 = 3. Also indi- 
cated in the figure are the masses of the objects as explained 
in the caption. We note that the star formation rates appear 
to be rather insensitive to the adopted value of r, indicating 
that for most of the simulated galaxies, the SFR* is rather 
con stant over a timescale of 0.1 Gyrs at high redshift (see 
also IWeinberg et al.l l|2003) for similar comparisons). There 
is some dispersion in the two estimates, mainly due to irreg- 
ularities in the star formation events in a given galaxy. In 
addition, some galaxies that have a non-zero SFR* as com- 
puted with r = 0.1 Gyrs, show no star formation on shorter 
timescales. These galaxies are also shown in fig. |21 but as- 
signed an arbitrary value of SFR*(t = 0.05) = 0.03 Mq/jt. 
For the low-mass objects with the lowest SFR* , there is an 
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SFR (t=10" yr) 

Figure 2. Comparison between the instantaneous star formation 
rates, in units of Mq/jt, computed for time dilated values of 
T = 0.1 Gyrs and t = 0.05 Gyrs for the catalog at z = 3. Galaxies 
with SFR* = 0.0 at T = 0.05 Gyrs are plotted at SFR* = 0.03. 
Symbols denote the mass range of each galaxy: 5 ■ 10* — 5 ■ IO^Mq 
(crosses), 5 ■ lO'' - IO^Mq (squares), M > IO^Mq (stars). The 
inset shows a blow-up of the highest SFR*-va,nge. 



apparent trend for the lower value of r to yield a lower value 
of SFR*, up to a a factor of a few. The larger value of r 
samples a longer period of star formation in the galaxy prior 
to the epoch under consideration and may therefore include 
periods of stronger star formation activity, resulting in a 
higher SFR*. 

Figure 13 shows the SFR* as a function of the mass of 
the objects for the catalogs at z = and z = 3. Galax- 
ies with non-zero SFR* axe plotted arbitrarily at SFR* = 
0.01 Mq/jt. In the former catalog, 572 objects have non- 
zero SFR* for r = 0.1 Gyrs, while in the latter these objects 
are 1470, but for clarity only 1000 randomly selected objects 
with non-zero SFR* are plotted. Also indicated in fig. |3 
are three ranges of the epoch of formation of the galaxies. 
For the catalog at z = 0, the dividing epochs are z = 0.9, 
and z — 1.1, roughly corresponding to half and the third of 
the age of the Universe. For the catalog at 2 = 3, we use 
z = 3.9 and z = 4.3 as the dividing epochs. For both red- 
shifts, the instantaneous star formation rate increases lin- 
early with mass although there is considerable dispersion 
especially at the lower redshift. For z — 0, this linearity 
holds only in the approximate mass range 10^" — 10^^i\'fQ, 
for lower masses there is a large scatter in the SFR*. For 
the higher redshift the linearity extends over two decades in 
mass. It is clear at both redshifts that the older objects tend 
to have lower SFR* than the younger objects and a slightly 
steeper mass dependence than linear, though the difference 
decreases with increasing redshift. 

Note also, that galaxies of a given mass tend to have up 
to 10 times higher SFR* at 2 = 3 than at 2 = and that 
the maximum SFR* in the catalogs decreases with redshift 
from about 100 Mq/jt at 2: = 3 to about 10-20 Mq/jt at 
2 = 0. It is interesting to note that these high SFR* val- 
ues at high redshift are consistent with the obse rvationally 
inferr ed SFR of Lyman-break galaxies at 2 = 3 llGiavaliscol 
|2003). The evolution of SFR* with redshift is very weak 
for the low-mass objects in fig. |3 as their SFR* tends to 
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Figure 3. Instantaneous star formation rate as a function of 
galaxy mass, computed with r = 0.1 Gyrs, for the catalogs at 
z = 3 (upper panel) and z = (lower panel). For clarity, only 
1000 objects with non-zero SFR* are plotted from the catalog 
at 2 = 3. The symbols indicate the epoch of formation of each 
galaxy. In the upper panel they denote: epoch of formation less 
than z = 3.9 (dots), epochs in the range z = 3.9 — 4.3 (triangles) 
and epochs greater than 4.3 (open circles). In the lower panel we 
have: less than z = 0.9 (dots), in the range z = 0.9—1.1 (triangles) 
and greater than 2 = 1.1 (open circles). In both panels galaxies 
with SFR* = 0.0 are plotted at SFR* = 0.01 Mg/yr. 



stay constant in time. The most dramatic difference in the 
cosmological evolution of the SFR* is seen in high-mass, 
old galaxies. The oldest (open circles) and most massive ob- 
jects have a SFR* up to two orders of magnitudes lower 
than their younger counterparts (dots). In fact, the star for- 
mation rate of massive, early-formed galaxies at 2 = is 
similar to that of the population of low-mass objects (few 
times 10^ Mq or less), almost all of which are late- formed 
galaxies. 

Note that the plots in fig. |3| show a number of galax- 
ies which are not forming stars at the time of observation, 
tz- In these objects star formation has already ceased or 
may be temporally stopped. These objects, although inac- 
tive, are quite numerous with most of them being of low- 
mass at both redshifts but being of higher mass as the red- 
shift decreases. At low redshift, most of the massive galaxies 
show no activity at all. These high-mass, inactive galaxies 
reside in the center of the most massive dark matter halos 
jAlimi fc Courtv^l2004^ and, as the evolution progresses, the 



GRB Host Galaxies 5 




log(10"yi- SFR'/M). 



(3) 



10^ 




10^° ^q11 
Mass (Mq) 






° .■ ' 




A . . 

AO. A 

OA A 

A 




o 






.■-AA 


z=0 r 



0.01 



0.1 



10 



100 



SFR (M„ yr ^) 



Figure 4. Upper panel: Efficiency parameter e, as a function of 
galaxy mass for the catalog at z = 0. Compare also with fig. 1^ 
where the different symbols are explained. Lower panel: e as a 
function of SFR* for galaxies in the catalog at z = 0. The dotted 
lines indicate a constant mass of 10'', IQl" and lO^M© (from left 
to right). 



depth of the potential wells prevents more and more gas from 
coohng. This hierarchical character of the galaxy formation 
process results in a strongly reduced star formation activ- 
ity at low redshift, and resembles that observed in elliptical 
galaxies in the center of galaxy clusters. Adopting a higher 
value of r would indeed cause some of these massive, early- 
formed galaxies to have a non-zero SFR* (recall fig. |5J as 
it would take into account the latest formed stellar popula- 
tions, in the case where there is still some cold gas available. 
This inactive population is important for our understanding 
of galaxy evolution and will be discussed further in Section 

m 

The amount of gas turned into stars at any given epoch, 
tz, depends on the amount of cold gas available in a gas 
region, hence the scaling of the SFR* with the mass of a 
galaxy. However, a low-mass galaxy can be more efficient in 
turning gas into stars than a high-mass galaxy. A measure 
of this efficiency is given by the specific star formation rate, 
the ratio of the SFR* to the mass of the galaxy. The specific 
star formation rate is usually normalized to a given mass, 
and we define a star formation efficiency parameter e by. 



log 
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Although e is a dimensionless parameter we will use it and 
the specific star formation rate, interchangeably and with 
the same meaning. In fig. 21 we plot e as a function of mass 
in the upper panel and as a function of SFR* in the lower 
panel, for galaxies in our catalog at z = 0. In the lower panel, 
galaxies of a given mass are located along diagonal lines, 
three of which are shown. It is apparent that the low mass 
galaxies have a much higher specific star formation rate than 
intermediate or high-mass galaxies and are therefore much 
more efficient in turning baryonic mass into stellar mass. In 
fact, the old and very high- mass objects, M > 10^^ -^0) s^nd 
some of the older intermediate mass objects, have specific 
rates that are up to two orders of magnitudes lower than the 
low-mass objects (A/ < 10^" -^0)- The results presented in 
the upper panel show a similar tre nd as seen in observat ional 
studies on blue compac t galaxies jGuzman et al Jll997f) and 
star- forming galaxies jpirez- 001^^1^^^1X1200311 . where 
the specific star formation rate is seen to increase with de- 
creasing galaxy mass for low-mass galaxies, M < 10^° Mq. 
Despite approximations and a number of uncertainties in- 
herent in numerical simulations, we find that our specific 
rates agree quite well with the observationally inferred val- 
ues. The l ow-mass objects are o bserved to have e in the range 
1.25 - 3 ||G uzman et al ]'1997), where we find the range to 
be about 1 — 2.5 in the simulations. 

The lower panel of fig. |1] reveals several interesting 
properties of the galaxy population. It is apparent that 
the galaxies tend to cluster in different regions in this fig- 
ure. Most of the objects in the intermediate mass range, 
have an efficiency in the range e = — 1, although there 
is a clear difference between the young and the old ob- 
jects, the younger ones in general being more efficient. The 
star formation activity of these intermediate mass objects, 
as defined by their star formation rate, spans the range 
SFR* = 0.1 — 20 Mq/jt. Another interesting feature seen 
in fig. 1^ are the two distinct groups of objects, one of high- 
mass with low efficiency (e<0), and the other of low-mass 
but with high efficiency (e > 1). The low-mass objects, that 
cluster around the diagonal line representing a galaxy mass 
of 10® Mq, have the highest specific star formation rate, but 
only an SFR* in the approximate range 0.1 — 1 Mq/jt. The 
most massive objects, although again spanning almost two 
orders of magnitude in their SFR* , distinguish themselves 
by having by far the lowest specific star formation rate. In 
general terms, we thus see that the low-mass objects are the 
most efficient in their star formation, by a factor of 10 over 
the intermediate class. The most massive ones are again a 
factor of at least 5 — 10 less efficient than the intermediate 
ones. 

We present in figurej^a statistical quantification of the 
star formation efficiency of the galaxy population at 2 = 0, 
where we have computed the number of objects in a bin of 
e. The dotted histogram in both panels represents all galax- 
ies in the catalog. This distribution peaks around e ~ 0.7, 
with considerable tails extending both above and below the 
peak. In the upper panel, the solid histogram shows the cor- 
responding distribution for all objects more massive than 
5 ■ 10^° Mq. The distribution is now cut off above e « 1.0, 
clearly demonstrating that it is the lower mass range of 
objects that have the highest star formation efficiency at 
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Figure 5. Distribution of galaxies per bin of e for galaxies in 
the catalog at 2; = 0. The dotted histograms, representing all 
galaxies in the catalog, are identical in both panels. The solid 
histogram in the upper panel is for galaxies more massive than 
5 • 10^" while the solid histogram in the lower panel is for 

objects with epoch of formation less than z = 0.9. All histograms 
are normalized to the number of objects in each distribution. 
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Figure 6. Efficiency parameter, e, of galaxies in the catalog at 
2 = 3 as a function of mass (upper panel), and as a function of 
SFR* (lower panel). For clarity, only 1000 objects are plotted in 
each panel. As in fig. |3] the dots denote formation epochs lower 
than z = 3.9, triangles denote epochs in the range z = 3.9 — 4.3, 
and circles denote formation epochs higher than 4.3. The dotted 
lines in the lower panel indicate a constant mass of 10®, 10^" and 
10"AfQ (from left to right). 



present. The solid histogram in the lower panel of fig.|S]dis- 
plays the corresponding distribution if only recently formed 
galaxies are included (epoch of formation less than z — 0.9). 
The main difference between the young population and the 
total is that the objects with low star formation efficiency 
have disappeared. Comparing the two panels, it is apparent 
that the low mass objects have the highest star formation 
efficiency (upper panel) and the lowest efficiency is generally 
found in the oldest objects (lower panel). In the following 
discussion we will denote the peak value of a given catalog by 
e, and refer to objects with e > e, as efficiently star-forming 
galaxies. Similarly, galaxies with e < e will be termed ineffi- 
cient. 



4 COSMOLOGICAL EVOLUTION 

The star formation activity of the galaxy population evolves 
strongly with redshift as seen in fig. |3 This evolution is also 
manifest in the cosmic evo lution of the star formation rate 
density (see lHopkins et all l|2001) for a compilation of obser- 
vational data) . Observations show a strongly decreasing star 



formation activity towards lower redshifts for z ^ 2, while 
its evolution towards high redshift is much less certain. The 
actual star formation rate in a given galaxy not only depends 
on the history of the galaxy, but also on merger events, en- 
vironment, etc. In this section we explore the evolution with 
redshift of the star formation properties discussed in the 
previous section. As the cosmic star formation activity de- 
creases for redshifts below z = 2, we consider mainly the 
catalogs of objects a.t z — and 2 = 3 for this exploration. 

In fig.|Hl we show e as a function of mass and star forma- 
tion rate for the galaxies in the catalog at z = 3. This should 
be compared to fig. |1| but note the different vertical scale. 
The efficiency is essentially constant over most of the mass 
range and the greatest dispersion is for objects with mass 
less than a few times 10^ Mq. The distribution of e is much 
narrower than at 2; = 0, and the mean value is also consider- 
ably higher, e « 1.8, an indication of the strong cosmological 
evolution of the efficiency of star form ation. This trend has 
also been seen in observa tional data iCuzman et al.lfl997l : 
iBrinchmann In the lower panel in fig. El we see 

that the clear separation of the galaxies into distinct groups 
or regions observed at 2 = 0, is nearly absent at 2 = 3. 
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Figure 7. Conditional moan of liaving a specific star formation 
rate at a given mass, {SFR/M\M), as a function of mass at differ- 
ent redshifts: z = (squares), 2 = 1 (diamonds), z = 2 (crosses), 
2 = 3 (triangles), 2 = 4 (stars), 2 = 5 (open stars); for clarity, er- 
ror bars representing dispersion around the mean are only plotted 
at 2 = and 2 = 3. 



The diagram does show the presence of a sub-population of 
low-mass objects with high efBciency but a moderate SFR* . 
This population is not clearly separated from the interme- 
diate mass objects at this redshift, but it continues to be 
present at z = with only slightly lower efBciency and ac- 
tivity. Comparing with fig. 21 we note another interesting re- 
sult, the appearance at low redshift of the high-mass objects 
with low efficiency and low SFR* . Note also that at both 
redshifts, the highest mass objects always have the highest 
SFR*. 

In fig. Q we show the conditional mean of having a spe- 
cific star formation rate at a given mass as a function of 
mass and its evolution with redshift. So as not to clutter 
the diagram too much we only show the dispersion around 
the mean for the catalogs at z = and z — 3. In general the 
dispersion is less at higher redshifts except for the low-mass 
objects as already noted in fig. 10] The largest difference in e 
between low and high-mass galaxies is seen at low redshift, 
where only the low-mass and some of the intermediate mass 
objects have high specific star formation rate. At high red- 
shift the entire galaxy population has a high specific rate. It 
is apparent from the figure that the specific star formation 
rate of the lowest mass objects decreases only slightly from 
z = 5 to z = 0, whereas for the intermediate to high-mass 
objects it decreases by two orders of magnitudes or more. 
It is clear that the low-mass objects are not only the most 
efficient ones at z = as already noted, they are efficient at 
all redshifts. At higher redshifts the difference between the 
efficiency of the different mass ranges diminishes, but the 
star formation efficiency of the low-mass objects is consis- 
tent with being constant out to a redshift of 5. 



5 DISCUSSION 

Numerical simulations provide us with catalogs of galaxy- 
like objects characterized by different properties such as 
mass, epoch of formation and star formation rate. Galaxy 
sub-populations in the catalogs are clearly distinct especially 
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Figure 8. Upper panel: Galaxy sub-populations at 2 = 1 as a 
function of mass and epoch of formation. Stars denote inactive 
galaxies (SFR* = 0), and squares active galaxies (SFR* > 0), 
with empty symbols for inefficient and filled ones for efficient 
galaxies (at this redshift the dividing efficiency is e = 1.3). Lower 
panel: Star formation efficiency, e, of the active population in the 
catalog at 2 = 1 as a function of SFR*. For clarity only 1000 ran- 
domly selected objects are plotted in both panels. As in previous 
figures, the dotted lines indicate constant masses of 10^, 10^0 and 
W^'^Mq (from left to right). 



at low redshift, e.g. the high-mass population with low star 
formation rate contrasting with the efficient low-mass pop- 
ulation. The former population shows a decreasing star for- 
mation activity with lower redshifts. As mentioned in section 
|3] a number of objects at low redshift have a zero value of 
SFR* when measured over a time r. To make the picture 
of our simulated galaxy population complete, we show in 
the upper panel of fig. |S] the entire population ai z — 1 
as a function of mass and epoch of formation. In this sec- 
tion we focus on the redshift 2 = 1, as it is close to the 
peak of the current GRB redshift distribution and most of 
the observed host galaxies have redshifts around unity. The 
catalog at 2 = shows the same trends as is easily veri- 
fied by comparing the lower panels in figs. |H] and 2] In fig- 
ure |H| stars refer to the inactive galaxy population that has 
no measurable SFR* at this redshift. Most of this popula- 
tion is either low-mass galaxies or high-mass galaxies. We 
note also that these objects are quite old. Among the active 
galaxy population (square symbols), galaxies with a non- 
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Figure 9. Number of galaxies per bin of mass at redshift of 2 = 3 
(upper panel) and z = 1 (lower panel). In each panel the solid 
histogram shows the distribution of the efficient population (with 
e > e = 1.8 at 2 = 3 and e > e = 1.3 at 2 = 1) and the dotted 
histogram shows the non-efficient population (with e < 1.8 at 
2 = 3 and e < 1.3 at 2 = 1). Both histograms are normalized on 
the whole population. 



zero SFR* (1239 objects) are separated in the upper panel 
of fig. IHl into an efficient population, including all galaxies 
with e > e — 1.3 (filled squares), and an inefficient popu- 
lation with e < e = 1.3 (open squares). The most efficient 
galaxies are low-mass and late-formed galaxies as already 
pointed out in the previous sections. The lower panel of fig. HI 
focuses on the active population in a e vs. SFR* diagram. 
This plot (as the one at z = in fig. 0J reveals again that 
high star formation rate galaxies are not necessarily galaxies 
with the highest specific SFR* . This kind of representation 
emphasizes the distinction between the star formation activ- 
ity {SFR* , along the horizontal axis) and the star formation 
efficiency (e, along the vertical axis), as two different mea- 
sures of the process of converting baryons into stars. 

The efficient and inefficient populations are statistically 
characterized at high and low redshifts in fig. 1^1 where we 
show their distribution per bin of mass. Note that both dis- 
tributions at each redshift are normalized on the whole pop- 
ulation, so adopting lower or higher values of e would alter 
the shape of both the efficient and inefficient populations. 
At 2 = 3, most of the objects with star formation efficiency 
above e = 1.8 (solid histogram) are low-mass objects (M < 



a few times 10 Mq). High- mass objects are not very nu- 
merous at this redshift (cf. fig. [SJ, and most of them are 
inefficient. At z = 1 most of the efficient objects are again 
predominantly of low-mass, although a considerable fraction 
of intermediate mass objects (approximately 10^*^ — IO^^Mq) 
also have high efficiency. The non-efficient objects at z = 1, 
are dominated by the mass range of few times 10^ Mq to 
almost lO'^^ Mq. At the lower redshift, essentially all of the 
low-mass objects are efficient in their star formation. The 
same can be said about the corresponding distribution at 
z — 0, almost all low-mass objects are efficient. To sum- 
marize, we note that efficient galaxies exist over almost the 
entire mass range at high redshift (see also fig. |7|l. At low 
redshift the efficiency becomes strongly mass dependent, as 
is also clear from the lower panel in fig.|3] Most importantly, 
the efficient population is dominated by a large number of 
low-mass objects, independent of the redshift. Considering 
the lower panels of figs.|H]and|H] we clearly see that the non- 
efficient objects have a larger spread in their SFR* than the 
efficient galaxies. The SFR* of the latter sub-population are 
clustered at moderate values. 

The star formation efficiency, e, of a given galaxy is ex- 
pected to be closely related to its color. Indeed, the instanta- 
neous SFR is dominated by the youngest stellar populations 
and thus is a measure of the luminosity in the optical/UV 
bands. The galaxy mass is on the other hand a better mea- 
sure of the luminosity in the infrared band, assuming some 
mass-to-light ratio, since the mass is the integration of all 
stellar populations produced over the lifetime of the galaxy. 
Low-mass galaxies of the simulated catalogs have the highest 
specific SFR* and are thus expected to be bluer than mas- 
sive galaxies; observations of compact galaxies show bluer 
colors a nd lower mass than high-mass disk-like starburst 
galaxies iGuzman et alJl^997^ . 

The first observational results on GR B host galaxies 
provid e important clues to their nature. iLe Floc'h et alJ 
i2003r showed that the hosts are f aint in the -RT band, typi- 
cal luminosities being about O.lL^. lBereer et al.l J2003t) have 
shown tha t the hosts are bluer than other sub-mm galaxies. 
ICharv eT^l. (_2002) considered 11 GRB host galaxies and 
derived extinction corrected SFR for seven of these hosts 
as well as mass, starburst age and internal extinction. Six 
of them have redshifts between 0.7 and 1.6, one is at a high 
redshift of z = 3.4. They derive UV and /3-slope SFRs in 
the range 0.2-70 Mq/jt for the hosts with redshifts around 
unity and more than 300 Mq/jt for the host with highest 
redshift. Using population synthesis spectral energy distri- 
butions they computed the mass of host galaxies and found 
them to be of rather low-mass, between 10** and 4 ■ 10^ Mq 
for the hosts with redshifts around unity. The host galaxy at 
z = 3.4 is a higher mass object, 3 • IO^^Mq. These mass vaf 



ues are also consistent with other estimates iSokolov et alJ 
200 li) . As a result, .Cha rv ot al. (2002) estimate the spe- 
cific star formation rate of these hosts, approximately in 
the range e = 2.5 — 3.9 in our units (see their fig. 2). Even 
with some uncertainties, these values are consistent with 
values inferred for star-bursting galaxies (see upper panel of 
fig . 16 in jPc rcz- Gonzalez et al. ( 2003) based on observations 
bv lCalzettil lll997l) ). Future observations will reveal more de- 
tails of the host galaxies and especially the connection these 
have with the observed population of blue galaxies, which 
could be irregulars, star-bursting, compact or Hll. 
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Based on these host galaxy observations and on the fact 
that we can identify galaxy sub-populations in the simula- 
tions, we define our simulated candidate GRB host galaxies 
as efficiently star-forming objects, with high specific star 
formation rate (filled squares in fig. (HJ A galaxy is consid- 
ered a host if its efficiency is higher than e, taken as the peak 
value of a given catalog. This value is rather arbitrary and is 
likely to be a lower limit, according to the observati onal esti- 
mate of the specific star formation rate in Chrlatensen et alJ 
i2004l) . The majority of these efficient galaxies in the sim- 
ulations are low-mass objects, but a modest peak in the 
mass distribution around a few times 10^° Mq is also seen 
in the lower panel of fig. |5| Such galaxies are expected to be 
arou nd 10 times less l umino us than L* galaxies, consistent 
with iLe Floc'h et al.l (|200.'^1 . Moreover, the simulated host 
galaxy candidates are moderately active, with SFR* span- 
ning two orders of magnitude around unity. Their specific 
star formation rates spans about one order of magnitude 
(fig. ISJ. Most of these galaxies are also low-mass objects, 
clustering around the line representing IO^Mq. This pic- 
ture seems consistent with observations: Although the host 
galaxies in ^harv ct al. ( 2002) show much higher specific 
star formation rates than our candidates, it varies by a fac- 
tor of 20 and the host masses are generally around 10^ -^Q- 
It is also interesting to note that one of their hosts with 
z = 3.4, has a similar e as the hosts around z — 1, but with 
a much higher star formation rate and a galaxy mass that 
is two orders of magnitude higher. The simulations show a 
similar result, a lO'"* Mq galaxy has a specific star forma- 
tion rate in the range e ~ 1.3 — 2.6 a.t z — 1 (see fig. IHJ, 
which brackets the values obtained for a lO'^" Mq galaxy at 
2 = 3 (fig.|nj. Recall also, that the SFR* increases at higher 
redshift for a given mass. 

The location of the candidate host galaxies in the 
e — SFR* parameter space, appears to indicate that most 
of them would have high e and mo dest to high SFR* . 
This is also what the observations of ICharv et alJ (l2002l) 
show. Moreover, a number of hosts are classified as compact 
galaxies, and such galaxies appear to have SFR between 
0.1 an d 14 Mq/jt as in the sample of iGuzman et alJ 
il997fl . Putting this together suggests that identifying host 
galaxies of GRBs as efficient galaxies provides us with first 
interpretations of the observations. Our results also point 
to a verifiable prediction regarding the observationally 
determined distribution of the specific star formation 
rates of GRB hosts. Observations of host galaxies could, 
however, be biased towards the most luminous objects and 
therefore the intermediate to high-mass range, with the 
more numerous and efficient, but much less massive, objects 
underrepresented. The solid histograms in fig. llUl show the 
efficiency distribution of the candidate host galaxies (i.e. 
objects with e > e) at high and low redshifts. Considering 
all efficient objects results in broad distributions (solid 
histograms). If we now restrict the hosts to be the most 
massive of the efficient objects (M > 10^'' Mq), we obtain 
strikingly different distributions. They are much narrower 
(dotted histograms) at both redshifts with the tail towards 
higher efficiency missing. First versions of observationally 
determined dis tribution s similar to the dotted ones in fig. 1101 
already exist iChristen scn 20021). The difference between 
the histograms, illustrate how our perception of the hosts 




log(10"yr SFR*/M) 

Figure 10. Upper panel: Number of host galaxies (e > e = 1.8) 
per bin of specific star formation rate (solid histogram) compared 
to the distribution of host galaxies with M > 10^*^ Mq (dotted 
histogram) at 2 = 3. Lower panel: same but at z = 1 with e > 
e = 1.3. Histograms arc normalized to the number of objects in 
each distribution. 



will change if we are able to detect these low-mass galaxies. 

The conclusion that host galaxies are, in majority, low- 
mass, late-formed galaxies with moderate star formation 
rates, derives from our hypothesis that candidate host galax- 
ies are the most efficient ones. We adopted this hypothesis 
as observations indicate that the hosts have high specific 
star formation rates and moderate SFR iChristensen et al] 
[2004). Our conclusions would be different if the hypothesis 
is that candidate hosts are galaxies with the highest star 
formation rates, as may be expected if the GRB formation 
rate follows the global SFR. In that case, the majority of the 
host galaxies would be high- mass, inefficient, early- formed 
galaxies (see Fig.lHJ. It may be argued that for our adopted 
hypothesis to hold, and if GRBs trace the global SFR, low- 
mass galaxies should contribute significantly to the global 
star formation rate. However, inspection of Fig. |H1 imme- 
diately shows that the matter is more complicated. If we 
select a sample of objects according to a particular criteria, 
we find that at z = 1, 50% of the total star formation rate 
is in galaxies with a SFR* < 9 Mq/yt; or in galaxies with 
e > 1.3; or in galaxies with a mass less than 5 ■ 10^° Mq. 
Similarly at z = 3, these thresholds become SFR* — 17 
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Mq /yr; or e = 1.9; or Af = 2 • IO^^Mq. Of course we have 
to keep in mind that these numerical values depend on the 
physics included in our simulations. 

It is clear from Fig. |S] that if a higher specific star for- 
mation threshold is used to select candidate host galaxies, 
fewer hosts will have high SFR*s. This would suggest that 
the GRB formation rate is not a tracer of the total star 
formation rate. The two rates follow each other (i) if the 
stellar mass function is independent of time and thus of the 
galaxy properties (such as mass, age, metal content, etc), 
and (ii) if the probability that a high-mass star gives birth 
to a GRB event is independent of time. Although a number 
of observational studies show t hat the stellar mass function 
depends on the g alaxy type iContini et alj [l995t iKroupal 
I2fl0ll : iLarsonlboort l , most workers assume a universal initial 
mass function independent of time. Moreover, the under- 
lying physics involved in the gamma-ray burst phenomena 
implies, among other things, that low-metalicity stars are fa- 
vored. Such low-metalicity environments should mainly exist 
in young or starbursting or interacting galaxies rather than 
in evolved systems even if these may be forming high-mass 
sta rs. Similar caveats h a ve be en dis cussed in other contexts 
by iRamirez-Ruiz et alJ ll2002l) and IChoudhurv fc Srianandl 
Observing host galaxies with properties such as a 
high specific star formation rate, could reinforce the fact 
that the probability of forming a GRB event from a high- 
mass star does depend on time and thus on the galaxy. 



6 CONCLUSIONS 

We have used numerical simulations of large scale structure 
formation to identify galaxy like objects and to follow their 
evolution. We concentrated on the mass, the epoch of forma- 
tion and the instantaneous and specific star formation rates 
of the objects. Our most important results are as follows; 

The whole galaxy population includes a significant num- 
ber of low-mass objects with M < W^°Mq. There is a 
strong cosmological evolution of the properties of the galaxy 
population between z — 3 and 2 = 0. We emphasize that 
the star formation rate of a galaxy should be considered a 
measure of its star formation activity, while the specific star 
formation rate is a measure of its star formation efficiency. 
These properties are nicely illustrated in a diagram of e vs. 
SFR* (fig. 0, where we find a number of clearly distinct 
sub-populations of galaxies. 

The low-mass population is further divided into sub- 
populations, some of which have a high specific star forma- 
tion rate. Based on these findings and on the first observa- 
tional results of GRB hosts, we identify candidate hosts in 
the simulations as galaxies with high specific SFR* , rather 
than objects with high star formation rate. The majority 
of these candidate host galaxies appear to be of low mass, 
with recent formation epochs and a moderate SFR* of a few 
Mq /year or lower. The properties of the candidate hosts so 
identified, are consistent with the trends observed in the 
GRB host galaxies. An implication of our conclusions, given 
the applicability of the physical input in the simulations used 
here, is that GRBs may not trace the cosmic star formation 
rate. 

Most of the observed GRB hosts to date have a redshift 
around unity. We show that the efficiency of the low-mass 



galaxies does not vary much with redshift. For high redshift, 
more and more galaxies are found to be efficient independent 
of mass, while the efficiency of intermediate to high-mass 
galaxies decreases strongly with decreasing redshift. 

As the low-mass objects are very faint, the observations 
of GRB host galaxies will initially be observationally biased 
towards the intermediate and high-mass objects. These will 
be easier to detect due to their higher luminosity. They will 
typically have an SFR in the range 1 — 40MQ/year at a 
redshift of unity. It may take the next generation space 
telescopes to reveal the properties of the low-mass sub- 
population of host galaxies. These are expected to outnum- 
ber the currently observable hosts by a considerable factor. 

Furthermore, the low-mass galaxy population is known 
to be important for the understanding of galaxy formation. 
GRB host galaxies may thus become useful tracers of this 
population or in fact the faint end of the galaxy luminos- 
ity function that may be difficult to detect otherwise. Thus 
GRB selected host galaxies may become an important link 
in the study of structure formation and evolution. Although 
GRB selected galaxies introduce its own 'selection effect', 
they are free from others afi'ecting surveys of various kinds, 
such as the magnitude limited observations of optical sur- 
veys and the sensitivity limitation of the sub-mm observa- 
tions. 
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